Modeling the Nd isotopic composition in the North Atlantic basin using an eddy-permitting model

International audienceBoundary Exchange (BE – exchange of elementsbetween continental margins and the open ocean) hasbeen emphasized as a key process in the oceanic cycle ofneodymium (Nd) (Lacan and Jeandel, 2005a). Here, we usea regional eddy-permitting resolution Ocean General CirculationModel (1/4) of the North Atlantic basin to simulatethe distribution of the Nd isotopic composition, consideringBE as the only source. Results show good agreementwith the data, confirming previous results obtained using thesame parameterization of the source in a coarse resolutionglobal model (Arsouze et al., 2007), and therefore the majorcontrol played by the BE processes in the Nd cycle onthe regional scale. We quantified the exchange rate of theBE, and found that the time needed for the continental marginsto significantly imprint the chemical composition of thesurrounding seawater (further referred as characteristic exchangetime) is of the order of 0.2 years. However, thetimescale of the BE may be subject to large variations as avery short exchange time (a few days) is needed to reproducethe highly negative values of surface waters in the LabradorSea, whereas a longer one (up to 0.5 years) is required tosimulate the radiogenic influence of basaltic margins and distinguishthe negative isotopic signatures of North AtlanticDeep Water from the more radiogenic southern origin watermasses. This likely represents geographical variations inerosion fluxes and the subsequent particle load onto the continentalmargins. Although the parameterization of the BEis the same in both configurations of the model, the characteristicexchange time in the eddy-permitting configuration issignificantly lower than the previous evaluations using a lowresolution configuration (6 months to 10 years), but howeverin agreement with the available seawater Nd isotope data.This results highlights the importance of the model dynamicsin simulating the BE proces

A modeling sensitivity study of the influence of the Atlantic meridional overturning circulation on neodymium isotopic composition at the Last Glacial Maximum

Using a simple parameterisation that resolves the first order global Nd isotopic composition (hereafter expressed as ε<sub>Nd</sub> in an Ocean Global Circulation Model, we have tested the impact of different circulation scenarios on the ε<sub>Nd</sub> in the Atlantic for the Last Glacial Maximum (LGM), relative to a modern control run. Three different LGM freshwater forcing experiments are performed to test for variability in the ε<sub>Nd</sub> oceanic distribution as a function of ocean circulation. Highly distinct representations of the ocean circulation are generated in the three simulations, which drive significant differences in ε<sub>Nd</sub>, particularly in deep waters of the western part of the basin. However, at the LGM, the Atlantic is more radiogenic than in the modern control run, particularly in the Labrador basin and in the Southern Ocean. A fourth experiment shows that changes in Nd sources and bathymetry drive a shift in the ε<sub>Nd</sub> signature of the basin that is sufficient to explain the changes in the ε<sub>Nd</sub> signature of the northern end-member (NADW or GNAIW glacial equivalent) in our LGM simulations. All three of our LGM circulation scenarios show good agreement with the existing intermediate depth ε<sub>Nd</sub> paleo-data. This study cannot indicate the likelihood of a given LGM oceanic circulation scenario, even if simulations with a prominent water mass of southern origin provide the most conclusive results. Instead, our modeling results highlight the need for more data from deep and bottom waters from western Atlantic, where the ε<sub>Nd</sub> change in the three LGM scenarios is the most important (up to 3 ε<sub>Nd</sub>. This would also aid more precise conclusions concerning the evolution of the northern end-member ε<sub>Nd</sub> signature, and thus the potential use of ε<sub>Nd</sub> as a tracer of past oceanic circulation

Global continental and marine detrital εNd: An updated compilation for use in understanding marine Nd cycling

Understanding the role of sediment-water interactions in the oceanic cycling of neodymium (Nd) isotopes is essential for its reliable use as a modern and palaeoceanographic tracer of ocean circulation. However, the exact processes that control Nd cycling in the ocean are poorly defined and require an up-to-date knowledge of the sources, sinks and transformation of this tracer to and within the ocean (e.g. as per the GEOTRACES core mission). We propose a considerable improvement of Nd-source identification by providing an extensive and up-to-date compilation of published terrestrial and marine sedimentary Nd isotopic measurements. From this database, we construct high resolution, gridded, global maps that characterise the Nd-isotopic signature of the continental margins and seafloor sediment. Here, we present the database, interpolation methods, and final data products. Consistent with the previous studies that inform our compilation, our global results show unradiogenic detrital Nd isotopic values (εNd ≈ -20) in the North Atlantic, εNd values of ≈ -12 to -7 in the Indian and Southern Ocean, and radiogenic values (εNd ≈ -3 to +4) in the Pacific. The new, high-resolution interpolation is useful for improving conceptual knowledge of Nd sources and sinks and enables the application of isotope-enabled ocean models to understand targeted Nd behaviour in the oceans. Such applications may include: examining the strength and distribution of a possible benthic flux required to reconcile global Nd budgets, establishing the potential use of Nd isotopes as a kinematic tracer of ocean circulation, and a general quantification of the non-conservative sedimentary processes that may contribute to marine Nd cycling

Resolving and parameterising the ocean mesoscale in earth system models

Purpose of Review. Assessment of the impact of ocean resolution in Earth System models on the mean state, variability, and future projections and discussion of prospects for improved parameterisations to represent the ocean mesoscale. Recent Findings. The majority of centres participating in CMIP6 employ ocean components with resolutions of about 1 degree in their full Earth Systemmodels (eddy-parameterising models). In contrast, there are alsomodels submitted toCMIP6 (both DECK and HighResMIP) that employ ocean components of approximately 1/4 degree and 1/10 degree (eddy-present and eddy-rich models). Evidence to date suggests that whether the ocean mesoscale is explicitly represented or parameterised affects not only the mean state of the ocean but also the climate variability and the future climate response, particularly in terms of the Atlantic meridional overturning circulation (AMOC) and the Southern Ocean. Recent developments in scale-aware parameterisations of the mesoscale are being developed and will be included in future Earth System models. Summary. Although the choice of ocean resolution in Earth System models will always be limited by computational considerations, for the foreseeable future, this choice is likely to affect projections of climate variability and change as well as other aspects of the Earth System. Future Earth System models will be able to choose increased ocean resolution and/or improved parameterisation of processes to capture physical processes with greater fidelity

The EC-Earth3 Earth system model for the Coupled Model Intercomparison Project 6

The Earth system model EC-Earth3 for contributions to CMIP6 is documented here, with its flexible coupling framework, major model configurations, a methodology for ensuring the simulations are comparable across different high-performance computing (HPC) systems, and with the physical performance of base configurations over the historical period. The variety of possible configurations and sub-models reflects the broad interests in the EC-Earth community. EC-Earth3 key performance metrics demonstrate physical behavior and biases well within the frame known from recent CMIP models. With improved physical and dynamic features, new Earth system model (ESM) components, community tools, and largely improved physical performance compared to the CMIP5 version, EC-Earth3 represents a clear step forward for the only European community ESM. We demonstrate here that EC-Earth3 is suited for a range of tasks in CMIP6 and beyond

Early to middle Eocene history of the Arctic Ocean from Nd-Sr isotopes in fossil fish debris, Lomonosov Ridge

Strontium and neodymium radiogenic isotope ratios in early to middle Eocene fossil fish debris (ichthyoliths) from Lomonosov Ridge (Integrated Ocean Drilling Program Expedition 302) help constrain water mass compositions in the Eocene Arctic Ocean between ∼55 and ∼45 Ma. The inferred paleodepositional setting was a shallow, offshore marine to marginal marine environment with limited connections to surrounding ocean basins. The new data demonstrate that sources of Nd and Sr in fish debris were distinct from each other, consistent with a salinity-stratified water column above Lomonosov Ridge in the Eocene. The 87Sr/86Sr values of ichthyoliths (0.7079–0.7087) are more radiogenic than Eocene seawater, requiring brackish to fresh water conditions in the environment where fish metabolized Sr. The 87Sr/86Sr variations probably record changes in the overall balance of river Sr flux to the Eocene Arctic Ocean between ∼55 and ∼45 Ma and are used here to reconstruct surface water salinity values. The ɛNd values of ichthyoliths vary between −5.7 and −7.8, compatible with periodic (or intermittent) supply of Nd to Eocene Arctic intermediate water (AIW) from adjacent seas. Although the Norwegian-Greenland Sea and North Atlantic Ocean were the most likely sources of Eocene AIW Nd, input from the Tethys Sea (via the Turgay Strait in early Eocene time) and the North Pacific Ocean (via a proto-Bering Strait) also contributed

The large-scale evolution of neodymium isotopic composition in the global modern and Holocene ocean revealed from seawater and archive data

Neodymium isotopic compositions (143Nd/144Nd or εNd) have been used as a tracer of water masses and lithogenic inputs to the ocean. To further evaluate the faithfulness of this tracer, we have updated a global seawater εNd database and combined it with hydrography parameters (temperature, salinity, nutrients and oxygen concentrations), carbon isotopic ratio and radiocarbon of dissolved inorganic carbon. Archive εNd data are also compiled for leachates, foraminiferal tests, deep-sea corals and fish teeth/debris from the Holocene period (< 10,000 years). At water depths ≥ 1500 m, property-property plots show clear correlations between seawater εNd and the other variables, suggesting that large-scale water mass mixing is a primary control of deepwater εNd distribution. At ≥ 200 m, basin-scale seawater T-S-εNd diagrams demonstrate the isotopic evolution of different water masses. Seawater and archive εNd values are compared using property-property plots and T-S-εNd diagrams. Archive values generally agree with corresponding seawater values although they tend to be at the upper limit in the Pacific. Both positive and negative offsets exist in the northern North Atlantic. Applying multiple regression analysis to deep (≥ 1500 m) seawater data, we established empirical equations that predict the main, large-scale, deepwater εNd trends from hydrography parameters. Large offsets from the predicted values are interpreted as a sign of significant local/regional influence. Dominant continental influence on seawater and archive εNd is observed mainly within 1000 km from the continents. Generally, seawater and archive εNd values form gradual latitudinal trend in the Atlantic and Pacific at depths ≥ 600 m, consistent with the idea that Nd isotopes help distinguish between northern/southern sourced water contributions at intermediate and deep water depths

Up to What Extent Can We Characterize Ocean Eddies Using Present-Day Gridded Altimetric Products?

The most common methodology used to detect and characterize mesoscale eddies in the global ocean is to analyze altimetry-based sea-level gridded products with an automatic eddy detection and tracking algorithm. However, a careful look at the location of altimetry tracks shows that their separation is often larger than the Rossby radius of deformation. This implies that gridded products based on the information obtained along track would potentially be unable to characterize the mesoscale variability and, in particular, the eddy field. In this study, we analyze up to what extent sea-level gridded products are able to characterize mesoscale eddies with a special focus on the North Atlantic Ocean and the Mediterranean Sea. In order to perform this task, we have generated synthetic sea level anomaly maps using along-track data extracted from realistic high-resolution ocean model simulations and applying an optimal interpolation procedure. Then, we have used an eddy detection and tracking algorithm to the gridded synthetic product and to the original model outputs and compared the characteristics of the resulting eddy fields. Our results suggest that gridded products largely underestimate the density of eddies, capturing only between 6% and 16% of the total number of eddies. The main reason is that the spatial resolution of the gridded products is not enough to capture the small-scale eddies that are the most abundant. Also, the unresolved structures are aliased into larger structures in the gridded products, so those products show an unrealistic number of large eddies with overestimated amplitudes

ORCHIDEE-ROUTING: revising the river routing scheme using a high-resolution hydrological database

International audienceThe river routing scheme (RRS) in the Organis-ing Carbon and Hydrology in Dynamic Ecosystems (OR-CHIDEE) land surface model is a valuable tool for closing the water cycle in a coupled environment and for validating the model performance. This study presents a revision of the RRS of the ORCHIDEE model that aims to benefit from the high-resolution topography provided by the Hydrological data and maps based on SHuttle Elevation Derivatives at multiple Scales (HydroSHEDS), which is processed to a resolution of approximately 1 km. Adapting a new algorithm to construct river networks, the new RRS in ORCHIDEE allows for the preservation of as much of the hydrological information from HydroSHEDS as the user requires. The evaluation focuses on 12 rivers of contrasting size and climate which contribute freshwater to the Mediterranean Sea. First, the numerical aspect of the new RRS is investigated, in order to identify the practical configuration offering the best trade-off between computational cost and simulation quality for ensuing validations. Second, the performance of the new scheme is evaluated against observations at both monthly and daily timescales. The new RRS satisfactorily captures the seasonal variability of river discharge, although important biases stem from the water budget simulated by the ORCHIDEE model. The results highlight that realistic streamflow simulations require accurate precipitation forcing data and a precise river catchment description over a wide range of scales, as permitted by the new RRS. Detailed analyses at the daily timescale show the promising performance of this high-resolution RRS with respect to replicating river flow variation at various frequencies. Furthermore, this RRS may also eventually be well adapted for further developments in the ORCHIDEE land surface model to assess anthropogenic impacts on river processes (e.g. damming for irrigation operation)